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Abstract 

By using simple mode coupling equations, we investigate the depo- 
larized light scattering spectra of two so-called "fragile" glassforming 
liquids, salol (phenylsalicylate) and CKN (Cao.4Ko.6(N03)i.4), mea- 
sured by Cummins and coworkers. Nonlinear integrodifferential equa- 
tions for the time evolution of the density-fluctuations autocorrelation 
functions are the basic input of the mode coupling theory. Restrict- 
ing ourselves to a small set of such equations, we fit the numerical 
solution to the experimental spectra. It leads to a good agreement 
between model and experiment, which allows us to determine how a 
real system explores the parameter space of the model, but it also 
leads to unrealistic effective vertices in a temperature range where 
the theory makes critical asymptotic predictions. We finally discuss 
the relevance and the range of validity of these universal asymptotic 
predictions when applied to experimental data on supercooled liquids. 
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1 Introduction 



Understanding the properties of the supercooled liquids is still a challeng- 
ing problem. This explains the considerable amount of work on the subject 
from the experimental as well as from the theoretical point of view. Of great 
interest is the fact that, whereas the static properties of a supercooled liquid 
vary weakly in the temperature range between the melting temperature T m 
and the so-called calorimetric glass transition temperature T g , the dynam- 
ical properties such as transport coefficients exhibit a strong temperature 
dependence. In particular, for 'fragile' liquids (according to Angell's classifi- 
cationB), the viscosity shows a strongly non arrhenian behavior : it increases 
slowly with decreasing temperature in the weakly supercooled domain (near 
T m ), but changes by many orders of magnitude when approaching T g . Below 
this point, the systems falls out of metastable equilibrium and becomes a non 
ergodic amorphous solid, the glass. Such an observation leads to assume the 
existence of two dynamical regimes crossing over at a certain temperature. 

From this experimental background, the mode coupling theory (MCT)i 
proposes a mechanism for the dynamics of supercooled liquids based on den- 
sity fluctuations : a closed set of integrodifferential equations is established 
within the formalism of the generalized Langevin equation by approximat- 
ing the memory function as polynomials of density correlators; a bifurcation 
scenario is then elaborated and describes the liquid-glass transition as one 
from an ergodic liquid state with zero Debye-Weller factor (DWF) to a non- 
ergodic ideal glassy state with non- zero DWF. A dynamical critical transition 
at which the singularity is expected to occur is then introduced at a tem- 
perature T c higher than T g . Due to its critical nature the MCT makes a 
number of universal predictions in the vicinity of T c . In particular, the in- 
dependence of the density correlators should vanish in the intermediate-time 
or -frequency domains where two relaxations, a and /3, overlap, the dynamics 
being described by a single function only depending on a unique parameter 
A. A square- root singularity of the DWF and a power- law divergence of the 
viscosity are predicted too. 
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To date most comparisons of the MCT to experimental data rely on these 
predicted universal features, although the range of validity of the approxi- 
mations leading to these predictions is not well determined (only recently, 
the necessary corrections to these asymptotic predictions have been investi- 
gated for hard spheres!). A different approach is to choose a reduced set of 
mode-coupling equations (then called "schematic equations") and to fit its 
full solution to experimental data with no direct reference to the expected 
critical properties.! Such a method, although very simplified, has the ad- 
vantage of taking explicitly into account the microscopic dynamics, which 
is impossible by other means, and of allowing an assessment of the range of 
validity of the universal predictions. A restricted study has already been pre- 
sented for the Raman spectra of several molecular glassformers (m-toluidine 
and other disubstitued benzenes, glycerol) in the frequency range from 150 
GHz to 6000 GHzJJ 

In this paper, we propose to apply this method to a wider frequency 
range! in order to fit the depolarized light scattering (DLS) spectra of salol 
(an aromatic molecular liquid) and CKN (a binary molten salt) between 
0.3 GHz and 5000 GHz, as measured by Cummins and coworkers. Many 
different experimental techniques have been applied to them, allowing tests 
of many aspects of the theory. A summary of these experimental results is 
proposed for comparison with our results in table [l] for salol and table |] for 
CKN. In the next section, we briefly describe the experimental support and 
some spectral features that our model has to take into account. In section 
[3], we review theoretical results needed for data analysis and discussion. Our 
model and data processing are presented in section |], whereas fitting results 
are exposed in section |5|. A discussion of the MCT predictions is given in 
section || and conclusions appear in section |7|. 

2 Experimental support and spectral trends 

We have analyzed the VH light scattering spectra of supercooled salol 
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and CKN obtained by Cummins and coworkers by combining both Raman 
and Brillouin scattering in the frequency range between 0.3 GHz and 5000 
GHz. The temperatures of experimentation cover the whole domain of the 
supercooled liquid state, from just above the melting point to the calorimetric 
glass transition. These spectra have been already extensively studied in the 
context of the universal predictions of the MCT (i'l for salol,l'0 for CKN) 
and the reader is referred to these articles for experimental details. 

Rather than the experimental scattered intensity, the susceptibility spec- 
tra are more suited for a comparison with MCT predictions. These are 
obtained by normalizing the experimental intensities by the Bose factor : 

„, , I exp (u) . , 1 . , 

n(uj) + l ex P(^r) - 1 

In the following, the term spectra will always refer to this frequency-depen- 
dent susceptibility. 

The spectra of salol present some general features shared by many aro- 
matic fragile supercooled liquids of the same kind like oterphenyl,00 m- 
toluidine,@ m-cresol and others. i These are : 

• a double high-frequency peak with principal frequency at about 3600 
GHz and a shoulder at 600 GHz, often referred to as 'Boson' peak 
and whose physical interpretation is not yet definitive.S&il This part 
of the spectra is almost temperature-independent and is assigned to 
microscopic excitations of the system. 

• a low frequency peak which maximum frequency decreases rapidly with 
temperature according to an non arrhenian law. It corresponds to the 
a-relaxation processes. 

• the intensity of the a-peak is a lot greater than that of the microscopic 
band. 

• the shape of the a-peak is characteristic of a stretched, rather than 
exponential, relaxation. 
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• a minimum occurs between these two structures with an enhancement 
of intensity at higher frequencies due to the so-called fast /3-relaxation : 
this region, where a- and /3-relaxations cross, is the domain of interest 
to test the universal predictions made by the MCT. 

The spectra of CKN, a fragile ionic glassformer, exhibit the same features 
except for two points : 

• the height of the a-peak is much lower, i.e. it has roughly the same 
intensity as the microscopic band. 

• no pronounced extra-structure is distinguishable in the microscopic re- 
gion. 

These differences may originate from the fact that both the translational 
and rotational molecular motions are detected, as suggested first for the o- 
terphenyl spectra.0 Indeed, the nature of the process that dominates VH 
spectra, interaction-induced or reorient at ion- induced scattering, has been in- 
vestigated for the different parts of the spectra of salol.il It was found that 
below the microscopic band the second mechanism is dominant, while a con- 
stant depolarization ratio of 0.75 was measured for all these liquids above 
600 GHz, supporting a dipole-induced dipole (DID) mechanism. Therefore, 
the observed differences could be simply related to a different weighting of 
these two contributions. This problem will be further discussed in section ||. 

3 Mode-coupling theory of the glass transi- 
tion 

As a starting point, the Mode Coupling Theory of the liquid-glass transi- 
tion assumes that the dynamics of a supercooled liquid is entirely governed 
by the density fluctuations and proposes a set of generalized kinetic equa- 
tions (d la Mori-Zwanzig) for the normalized autocorrelation functions 4> q (t) 
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of the density fluctuations modes at wave vector q, p q (t)$ These equations 
take the form : 

4> q {t) + n 2 4> q {t) + /* M q {t - 1') 4> q {t') dt' = o, (3.1) 



where the memory kernel is split into an instantaneous damping term 7 g 5(t) 
and a mode-coupling term m q (t) : 

M q (t) = n 2 g [ lq S(t)+m q (t)}. (3.2) 

The core of the mode coupling approximation in the "idealized case" only 
considered here is to express m q {t) as a polynomial function F q of the g (i)'s, 
so that 



(t) = F q (V,Mt)) = i: E V qii ..., qn <j> qi {t) ...<t> qn {t), (3.3) 



n>l qi,...,q„ 



where the components of the vector V, the vertices V qi} ,„, qn , are positive con- 
tinuous functions of the static structure factor S(q). It should be noticed 
that the linear terms in m q (t) are introduced in an ad hoc fashion to guar- 
antee that models with only a very reduced set of equations reproduce the 
stretched behavior of the a- relaxation. Some justifications of the presence of 
these terms have been proposed in the context of the non-linear fluctuating 
hydrodynamics!^! The mode-coupling equations are finally written as : 



ct> q {t) + lq Q 2 q 4> q (t) + Q 2 <f) q (t) + tt 2 / m q (t - t') ct> q {t') dt' = 0, (3.4) 







with m q (t) given by equation |3~3 . 



This set of equations can lead to an ergodic-to-non ergodic transition. 
Indeed, the Debye-Weller factors (DWF) defined by 

f q = lim <f> q (t) (3.5) 

t — » + oo 

are solutions of the set of equations : 

TZTT = F q (y,fk). (3-6) 
1 Jg 
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The only real solution in the vicinity of V = is the trivial one (for all 
g's, f q = 0) corresponding to an ergodic liquid state. A bifurcation to a non 
zero solution corresponding to a non ergodic "ideal glassy state" occurs when 
a critical hypersurface, the "ideal glass transition hyp er surf ace" , is crossed. 



Two type of transitions exists : the A-type with the DWF varying continu- 
ously through the transition hypersurface and the B-type with discontinuous 
change of the DWF. The second one is the only potentially relevant in the 
context of the liquid-glass transition. At the thermodynamic level, intensive 
critical parameters (critical temperature, density and so on), at which the 
transition is supposed to occur, are defined. 

In the vicinity of the transition hypersurface, it is possible to make various 
predictions on the critical behavior of the solutions of the mode coupling 
equations, the most important being the "reduction theorem". Expanding 
to lowest order in V — V c , one obtains 



where the leading contribution to G q , in the vicinity of the plateau of <j> q , is 
proportional to a g-independent function G, which is solution of the equation 
(called the scaling equation) 



It is defined by a vector V c such that the Jacobian matrix of equations [T6] 
is singular, so that : 




<f> q = r q + {i-r q ) 2 G q {t), 



(3.9) 



a + A G{tf - ^- f G(t - t') G(t') dt' = 0, 



dt Jo 



(3.10) 



where a, the separation parameter, is an 0(V — V c ) given by 




(3.11) 
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and is positive in the ideal glassy state, negative in the liquid, and zero on 
the transition hypersurface. The exponent parameter A is defined by 

X = \T.^ (w^f {YCJn ) (l-^a-^e*^, (3-12) 

q,k,p \ J k J p ) 

with [e q ] and [e q ] eigenvectors of the matrix |3l| l 

Eqs. p.9| and |3.10| and some of their consequences have been extensively 
used to test the validity of MCT by comparison with experimental data. 
In particular, an approximate expression around the susceptibility minimum 
has been proposed : 

Xmin (if ^ \ a i f^min^ 



a + b V Vcu m ; n / V a; 



with a and 6 such that 



r(i-a)* _ r(i+ff 
A " r(i-2a) " r(i + 2 6)- (3 ' 14) 

The relevance of such studies will be discussed in section || 



4 Model and data processing 

In our calculations, we intend to take into account the properties of the 
spectra in the THz domain such as the microscopic peak and the inelastic 
boson peak. This is not possible with studies by means of the scaling equa- 
tion, eq. |3.10| or eq. |3.13| , because they are restricted to the region around 
the minimum of the susceptibility spectrum, i.e. to the GHz domain. Thus, 
we choose to study in full detail a schematic model, i.e. a restricted set 



of mode-coupling equations like eq. [3]4| with a simple form of the memory 
functions, eq. |3~3|. 



The basis of our study is the well-studied F 12 -modeli for a single correlator 
i>o defined by 

Mt) + r o Mt) + fi o Mt) + fi o Cm{t - Mt') dt' = o, (4.1) 
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m (t) =v 1 Mt)+V2<t>l(t). (4.2) 

We add a second correlator 0i in order to reproduce the double-hump shape 
of the peak in the THz domain and the extra-intensity of the a-peak. It is 
solution of 

4>! (t) + Fx fii 0i (t) + n\ 0! (t) + Ql fmi (t - t') 0i (*') dt' = 0. (4.3) 

J 

It remains to express the memory kernel mi. For this purpose, we take the 
first correlator 0o accounting in an effective way for all modes that are re- 
sponsible for the underlying critical transition (typically, these are thought to 
be the density fluctuations at wave vectors corresponding to the maximum of 
the static structure factorial). The second correlator X describes additional 
degrees of freedom, that are important in depolarized light scattering pro- 
cesses but whose slow time-dependence is dominated by that of <fto (so that 
the liquid-glass transition occurs systematically for the two correlators at the 
same time). A simple choice is then 

mi(t) = rm (t). (4.4) 

In addition we impose that the characteristic pulsations (flo, the damp- 
ing coefficients (ro, Ti) and the coupling parameter r are temperature- 
independent. As a consequence, the whole temperature dependence is con- 
tained in the vertices V\ and v 2 - 

The above model has a B-type transition, and the critical hypersurface 
reduces to a line in the 2-dimensional parameter space (v±, V2) '■ 



v\ = 2 sjv c 2 - v c 2 with 1 < v c 2 < 4 (4.5) 
The exponent parameter is : 

A = (4.6) 

Other models with two correlators have already been proposed in a similar 
context : the Sjogren model for diffusion of impurities in a glass-forming 
liquid@ and the Bosse-Krieger model for coupled mass-density and charge- 
density correlators in a supercooled molten salt.il 
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It should be noted that, because of the extreme simplicity of the model, 
the real nature of the correlators is somewhat obscured. 0o and <pi are not 
necessarily restricted to be density correlators. In particular, the model could 
account for molecular orientional fluctuations or any other degrees of freedom 
that are not included in standard MCT. Accordingly, the parameters that we 
extract from the fits do not have a precise microscopic content and should 
be considered as a minimal set of phenomenological quantities that allows to 
describe the global evolution of the systems in the context of simple mode- 
coupling models. 

To compare the predictions of the model to the experimental data, it is 
also necessary to make some assumptions about the mechanism that gives rise 
to the light scattering spectra. When the dipole-induced dipole (DID) mech- 
anism dominates, as often assumedJIS the Stephen factorization schemeii 
(which is compatible with the factorization approximations underlying the 
MCT) leads to an expression for the DLS spectrum ~x"{oj) which involves 
the Fourier transform of bilinear products of the correlators. A simple de- 
scription of the spectrum for the case of our two-correlator model is similarly 
given by : 

X» = u A Im {FT [Mt) 2 + lMt) 2 } } (4-7) 

with FT denoting Fourier transform and Im the imaginary part. A is an 
amplitude factor that we take dependent of temperature. 7 is an ad hoc 
parameter which specifies the relative weight of the two contributions to the 
DLS spectrum due to 0o and to <p± and which we take as temperature inde- 
pendent. However, both A and 7 are substance-dependent. This expression 
of the scattered intensity corresponds to the assumption that 0o and 0i are 



both density correlators. Nevertheless, as pointed in ref. |22], fluctuations 
of the molecular orientations are a major contribution to the DLS spectra. 
In order to test the dependence of our results in the chosen expression of 
the scattered intensity, we have also performed tests with other choices for 
the expression of the susceptibility x"^); like a linear combination of the 
correlators or 'mixed' models, quadratic in O and linear in The fits 
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were not as good as with eq. |4.7| in the microscopic part of the spectra, but 
they were equivalent in the region of the minimum of susceptibility and of 
the a-peak. The values of the vertices and of other parameters of the model 
were only very weakly affected by the modification. It seems thus that the 
question of the nature of the light scattering mechanisms is not determining 
for our study. It remains out of our scope. 
The calculations consist of three steps : 

• 1) the time evolution of the correlators 0o and is computed from 
equations [4.1| to by means of an algorithm derived from the one 
described in reference^ and provided to us by M.Fuchs and W.Gotze. 



• 2) The Fourier transform in equation [4.7| is performed via the Filon 
algorithmic to obtain theoretical susceptibilities x"{ u )- 

• 3) The preceding calculations are included in a loop using a modi- 
fied least-square Powell algorithmil where the T-dependent parameters 
A,vi and t>2 and the T-independent parameters Qq, Q±, Tq, Ti, r and 
7 are adjusted to provide the best fit between model and experiment. 

Analogous calculations have been performed for different choices of the 
polynomial approximation to the memory function mo to fit Raman scat- 
tering spectra of m-toluidine.0 Good and stable fits were obtained with the 



F 12 -model associated with equation [17| for x"^)^ so we restrict mostly our- 



selves to this case. We have used different functional forms of the memory 
functions to fit the salol spectra, but the results are still preliminary and will 
be only briefly mentioned.^ 



5 Numerical results 

A major difficulty with the MCT in its formulation only in terms of 
density fluctuation modes is that it predicts the existence of a sharp ergodic- 
to-non ergodic transition that is not observed experimentally. This transition 
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which is assumed to be at a temperature T c higher than T g is avoided because 
of so-called activated processes that are not taken into account by the theory 
and that restore ergodicity below T c . For this reason, the low-temperature 
features predicted below T c are always obscured in the supercooled liquid 
state and the rare characterizations of theseB do not come without any 
polemics.ii It is thus necessary to impose some temperature limitations to 
the fitted data by ignoring the coldest spectra or imposing a low-frequency 
cut-off for the lowest temperatures in order to free our calculations from this 
difficulty without loosing any microscopic information. It is clear that such 
constraints will affect the fitting results. This point is discussed later on. 

Typical fits to experimental data are shown in figure |TJ at different tem- 
peratures for both liquids. With the chosen cut-offs (typically about 100 
GHz for temperatures below the best fit to T c . See fig. |l[a), the schematic 
spectra are in very good agreement with the data. In particular, at high 
temperature, the whole spectra are reproduced by the model including well 
defined a and microscopic peaks. Thus this set of equations in the mode- 
coupling formalism ensures the continuity of the dynamical processes from 
the transient vibrational regime of the microscopic band to the regime of 
structural relaxation. A drastic difference of behavior exists between high 
and low temperatures as far as a choice of frequency cut-off is concerned. 
Out of the restricted fitting frequency range, a reasonable extrapolation to 
low frequencies is possible for high temperatures, whereas it is impossible 
at low temperatures. This is due to the limitations of the theory discussed 
above. In the first case (high T), the spectra can be described unambiguously 
by the model, but, in the second case (low T), the theory does not describe 
the a-relaxation at all and the fits force the low-temperature predictions of 
MCT to occur below the frequency cut-off, giving unstable and unphysical 
results. 

The calculated temperature-independent parameters of the model are pre- 
sented in table |3|. The choice of frequency cut-off has little influence on these 
values. Error bars of about 15 % are determined from different choices. The 
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values obtained for the frequencies are realistic in regard to the position of 
the microscopic bands. 

The coupling parameter r measures the strength of the coupling of the 
correlator <pi to the driving correlator O - From the fits to both salol and 
CKN data, it is found greater than one. This supports the assumption that 
processes can be isolated that drive the overall dynamics of the system. 

An interesting result is that the relative weight 7 of the second correlator 
contribution to the DLS spectra is greater than one and greater for salol than 



for CKN. This is illustrated in figure where the two terms of equation 4.7 



are considered separately. It is clear from the figure that the main contribu- 
tion to the spectra is due to the second correlator except for the highest 
frequencies. Therefore the first correlator seems to drive the overall slowing 
down of the system, whereas the second one provides the experimental probe 
of the dynamics. It is an interesting consequence of the fitting procedure 
chosen with eq. \i.7\ to allow such a classification. Moreover, we may assess 
the contribution of the second correlator depending on the chemical nature 
of the liquid (aromatic or ionic), at least with our specific F12 model. Depart- 
ing from the original MCT scenario that deals only with density fluctuations 
and since our results are independent of the expression of the scattered in- 
tensity as a function of the correlators (see section ^), thus of their precise 
physical content, one could interpret this result in a way compatible with ref. 
f22| : orientational fluctuations represented by <pi produce the main part of 
the depolarized light scattering spectra and are strongly coupled to density 
fluctuations represented by O which drive the dynamics of the supercooled 
liquid above T c . It would be necessary to support this assumption to demon- 
strate that the dynamics of the orientational and density fluctuations obey 
similar mode-coupling equations, which remains to be done. 

We have also let r and 7 vary with temperature. However, it does not 
produce better fits despite the considerable amount of new adjustable pa- 
rameters and it was abandoned. 

The most important output of our calculations are the effective-temp er- 
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ature dependent vertices V\ and v% (figures [3] and f|). Indeed, the knowledge 
of these parameters gives access to the way a realistic system explores the 
parameters space of the model (see figure f|). It thus allows to discuss the 
applicability of the asymptotic predictions of the theory in the vicinity of T c . 
See below. 

For determining the vertices, the problem of the low-frequency cut-off 
must be addressed, again because of the limitations of the theory discussed 
above. At high temperatures, the vertices are determined unambiguously, 
irrespective of the cut-off : whether one restricts the fitting range to the mi- 
croscopic peak and its "fast-/3-relaxation" wing (i.e. a cut-off at about 100 
GHz) or includes the region of the susceptibility minimum, and considers 
or not the a-peak has no influence on the results. On the contrary, at low 
temperatures where the theory cannot apply to the full spectrum, the fitting 
procedure gives results that are cut-off- dependent and unstable. This pro- 
duces unphysical discontinuities in the behavior of v\ and v<i when crossing 
the dynamical-transition line. With no frequency cut-off (other than the one 
due to the experimental frequency window), the limitations of the theory 
appear in a rather surprising way on the fitted values of the vertices. At 
the lowest temperatures, the representative curves of salol and CKN seem to 
follow asymptotically the transition line as the temperature is lowered (see 
figure |]). Implications of this observation are discussed in the next section. 

Note finally that the choice for the expression of the DLS susceptibility 
as a function of the correlators 0o an d 4>i has no influence on the values of 
the vertices that we obtain. Linear or quadratic expressions produce undis- 
tinguishable results. This supports our analysis and conclusions, which are 
independent of the detailed description of the light scattering process. 

6 Critical properties 

Knowing how a real system explores the parameter space of the model 
can give some insight into the applicability of the MCT universal predictions 
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to experimental data. This is the purpose of our schematic model to provide 
such a knowledge. 

As discussed in the preceding section, the "idealized" MCT used here fails 
in describing the dynamics of the strongly supercooled liquids, i. e. below the 
critical temperature T c introduced by the theory itself. It is thus necessary 
in general to impose temperature and/or frequency limitations to the experi- 
mental data under study in order to compare with the theoretical predictions. 
Fortunately, this has no influence on the high-temperature results. For this 
reason, we consider in the following only the results obtained without explicit 
frequency cut-off. Doing so, the cut-off dependence of our results, a possible 
objection to mode-coupling tests,!! is avoided and the comparison with other 
tests of the theory is made easier. The penalty is of course that one must 
impose a low-temperature cut-off. 

The adjusted effective vertices exhibit two regimes with temperature. 
Above and slightly below the melting temperature, they vary linearly (v% is 
almost constant), whereas, at lower temperature, they tend to follow asymp- 
totically the transition line without crossing it. As the static structure fac- 
tor evolves only slowly with temperature in the supercooled liquid state, the 
rapid departure from the quasi-linear behavior appears to signal a breakdown 
of the usual MCT hypothesis that the vertices are continuous functions of 
the static correlations and vary moderately with temperature. If the criti- 
cal properties are determined by extrapolation from the quasi-linear domain 
only (table 01), the critical temperatures T c for salol and CKN are found in 
good agreement with the previous determinations (compare table Hj and ta- 
bles [I] and HD. It must be noted that the corresponding critical temperatures 
fall in the low-temperature but still ergodic regime. An exponent parameter 
Xht can be determined from this high-temperature behavior. It is in good 
agreement with previous determinations for salol but is significantly smaller 
for CKN (see tables). 

Other extrapolations are possible by including the nonlinear regime. Once 
one includes this regime, the critical temperatures and the exponent param- 
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eters depend on the low-temperature cut-off, and, if one wants to include 
more data, on low-frequency cut-offs as well. It is found that the smaller the 
low-temperature cut-off, the smaller the extrapolated critical temperature, 
and the smaller the exponent parameter A, as can be seen in figure |j. In 
an extreme situation, without any frequency cut-off (and with a wider ex- 
perimental frequency window), because the fitting procedure forces T c to be 
smaller than all the considered temperatures, it appears that no transition 
could be defined at all, as the vertices follow asymptotically the transition 
line without crossing it. Unphysical values could then be obtained, like v\ 
becoming negative. 

Of particular interest is the exponent parameter extrapolated from the ef- 
fective vertices at the critical temperature determined from the linear regime. 
It is referred to as \ vert in table |[ Whereas one would expect X vert to be 
compatible with the previous determinations of A based on the critical pre- 
dictions of the theory (tables [I] and |2|), it is found significantly smaller for 
both salol and CKN. This discrepancy should then affect strongly the quality 
of a master equation fit and limit its range of validity. 

We have also proceeded with different schematic models, a Fi 3 -based one 
with 

m = v 1 (p + v 3 (pi m 1 = rm (6.1) 
and a Sjogren-like onell with 

m = vi 4> + v 2 4>l m x = r <p <j) X . (6.2) 

Results of the fitting procedure are in total agreement with the above pic- 
ture M 

It is interesting to relate our results to the numerical studies on binary 
mixtures of Lennard- Jones atomdH and soft spheres.il For both systems, 
the static structure factors have been used to determine the vertices of the 
memory kernels in the mode-coupling equations. Then, the critical properties 
and the long-time dynamics obtained from the equations have been compared 
to the molecular dynamics simulation results. It is systematically found 
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that this mode-coupling approach overestimates the ability of the system to 
freeze in an " ideal glassy state" . The physical interpretation of this effect is 
clear : restricting the relevant dynamical modes to the density fluctuations 
excludes possible alternative decay channels for relaxations and the freezing 
occurs more rapidly. In the case of the Lennard- Jones system, the exponent 
parameter calculated from the mode-coupling model (A = 0.708; for soft 
spheres it is A = 0.73) is found smaller than that determined from a fit 
to the scaling equation (A = 0.78 ± 0.02). For both systems, the diffusion 
constant seems to vanish (at an extrapolated temperature) with a power law 
corresponding to A ~ 0.60, whereas the a-relaxation times seem to diverge 
with power laws corresponding to A ~ 0.79 for the Lennard- Jones system and 
A ~ 0.88 for the soft spheres. This is in contradiction with the requirement 
that a given system be characterized by a single exponent parameter A. 

These discrepancies on the value of A are surprinsingly quantitatively 
very similar to those found from our schematic calculations, although rig- 
orous and complete sets of mode-coupling equations (and not an arbitrary 
phenomenological model) are considered. Consequently, it seems reasonable 
to us to assume that a schematic study of the molecular dynamics simulation 
results would reveal effective behaviors similar to those found for salol and 
CKN. 

On the basis of the preceding discussion, we suggest a possible descrip- 
tion of the approach of the calorimetric glass transition in term of effective 
temperature-dependent mode-coupling parameters. At high enough temper- 
atures, the dynamics of a liquid are appropriately described by the ideal 
mode-coupling theory (in the sense that the dynamics are dominated by 
processes that couple the way the density correlators do in the MCT with 
vertices that are smooth functions of the static correlations in the system). 
When the temperature is lowered, the system seems to evolve toward the 
ergodic-to-non ergodic transition predicted by the theory. But, at a temper- 
ature notably higher than the extrapolated T c , the dynamics depart from the 
predictions of the theory : it is better expressed in term of effective vertices 
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that are not directly related to the statics and behave such that the transition 
is avoided. Finally, at low temperature, a qualitative change in the nature of 
the relaxation processes reveals a possible underlying dynamical transition. 
Quantitatively, this picture exhibits two properties : 

• the ideal mode-coupling approach leads to an overestimated ability of 
the liquids to freeze. 

• incompatible values of the exponent parameters are determined : com- 
pared to the one calculated from the statics, the fitting of the scaling 
equation overestimates A, whereas the effective behavior leads to an 
under est imat ion . 

If this picture works well for CKN, this last point is not clear in the case 
of salol, where our schematic fits provide equal exponent parameters from the 
high temperature extrapolation and the scaling equation. It is thought to be 
a peculiarity of this system, when approached with our elementary model. 
In CKN, there is a relative equivalence between the contributions of the two 
correlators to the scattered intensity (the proportions are about 30%- 70%), 
whereas in salol the second correlator represents about 95 % of the scattered 
intensity, as a consequence of the height of the a-peak much larger than the 
microscopic peak. Thus our calculations for salol reduce in reality in the 
GHz domain to a fit of this single correlator, that is in fact equivalent to a 
fit of the scaling equation in this domain (see figure |2]). Only more complex 
models with more correlators could reveal the proposed picture. 

The nature of the processes responsible for the transition to an effective 
mode-coupling behavior is of course of great interest. It is not clear that dis- 
crete processes as correlated jumps and hopping are good candidates, because 
in simulations they are generally characterized only well below the extrapo- 
lated critical temperature of MCT (for Lennard- Jones mixture,@ T c = 0.922 
is predicted from statics, whereas the secondary peak of the self part of the 
van Hove correlation function generally associated to these processes appears 
below 0.5). 
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If the picture described here is valid, it introduces some questions about 
the relevance of the predictions of MCT in the vicinity of the dynamical 
critical transition, that we found to be widely avoided. It is thus necessary 
to know what is the maximum distance to the critical hypersurface at which 
the universal predictions of the theory are effective : schematic models can 
provide some answers to this determining problem. 

First consider the elementary F 12 -model for one density correlator defined 
by equations (11] and Q.2\ alone. Following Gotze and Sjogren,^ we choose 



Vl = vl(l-e) V2 = v c 2 (l-e) (6.3) 

such that A = l/x/uf equals 0.7 on the transition line, e is thus a kind of 
'distance' from the transition line. After numerical integration and Fourier 
transform, the susceptibility spectra are fitted in the minimum region with 
the scaling equation (equation |3.10| ) in order to determine an effective ex- 



ponent parameter A e ^ as a function of positive e (figure ||). This procedure 
reproduces the usual way of treating experimental results beyond the scope 
of MCT. It is clear that for e greater than 0.005, X e ff departs notably from its 
value on the critical line and is never constant. It corresponds thus to a max- 
imum distance from the transition line for the critical properties predicted 
by the theory to be effective. 

In order to take into account this e dependence, a vertices dependent 
exponent parameter can be defined following Gotze and Sj6gren.il In the 
case of the i^-model and of our schematic model, which have exactly the 
same critical properties, this reads : 

hoc = v 2 (1 - /) 3 with / solution of v 1 + 2v 2 f = M 1 „, 9 (6.4) 

(i - fr 

valid only if 

vi > 3 vl /3 - 2 v 2 (6.5) 

We see that, in the F 12 -model, the introduction of local scaling equation 
parameters extends the domain of numerical compatibility of this equation 
to e = 0.015, but it is not enough : effective parameters are determined out 
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of this range and even in the domain where the condition |67| is not obeyed. 
These results have thus no significance in term of the MCT and arise only 
because of the numerical flexibility of the fitting procedure. 

The same problems should appear with the experimental spectra. Indeed 
the trajectories followed by salol and CKN in the parameters space of our 
model are not nearer, in term of a distance defined as e is, from the transi- 
tion line than are the points considered above. Moreover some experimental 



points do not obey the constraint |6T5 . 

We have performed scaling equation calculations as above for the suscep- 
tibility spectra of salol between 263 and 333 K. It is clear from figure |6| that 
the values of \ e jj vary notably with temperature and can be determined even 
when p.5| is not valid. They are compatible with the previous determina- 
tions by other authors (table P but neither with X vert (table |]) nor with Xi oc 
coming from our schematic calculations.^! This is representative of the inef- 
fectivity of the reduction theorem, as is the existence of two distinct minima 
for each correlator at all temperatures (figure Two possible reasons can 
be proposed for this, depending on the considered temperature domain. At 
high temperature, where the description of the dynamics by MCT is correct 
and alternative processes are not necessary, the supercooled liquid is too far 
from its transition point and thus the asymptotic predictions of the theory 
valid in the vicinity of the transition do not yet hold, whereas, at low temper- 
ature, processes that are not tractable within the mode-coupling formalism 
are probably at the origin of these discrepancies. 

The point is thus that our simple schematic calculations indicate that the 
leading order critical predictions of the MCT are experimentally unobserv- 
able. The conclusion would have been the same even if there were no need 
of additional processes to describe the dynamics of the supercooled liquid in 
the transition region; it is found necessary the system to be very near its 
transition point, where error bars of mosts experimental techniques impose 
differences of several degrees to reveal significant evolutions. This necessity 
has been pointed out recently within the investigation of the next-to-leading 
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order corrections to the critical predictions of MCT,S which are not for the 
moment available for experimental data like DLS spectra. 

7 Conclusion 

We have analyzed depolarized light scattering spectra of two fragile glass- 
forming liquids, salol and CKN, by means of schematic mode-coupling equa- 
tions (i.e. a reduced closed set of non linear integrodifferential equations for 
correlation functions). Such simple models are interesting because they make 
possible tests of the MCT that are not limited to its asymptotic critical pre- 
dictions. Moreover, they capture the essential structure of the mode-coupling 
equations with no reference to the details of the physical processes at work. 
It is thus possible to take into account in an effective manner processes that 
have not been included so far in the MCT. It is the case for instance of the 
orientational fluctuations. 

We have found that the two-correlator model developed here gives a very 
good account of the experimental data. In particular, the so-called boson 
peak offers no extra-difficulty within our fitting procedure, whatever func- 
tional form is chosen for the memory kernel.Bil 

The evolution with temperature of the relaxational properties of the sys- 
tem can be rationalized in term of effective temperature-dependent mode- 
coupling vertices. It is compatible under certain assumptions with observa- 
tions previously made in molecular dynamics simulations of simple liquids. 
Two temperature regimes can be introduced and interpreted as follow. At 
high temperature, the dynamics are well described by mode-coupling equa- 
tions, in which the vertices are determined from the static correlations of 
the supercooled liquid. This regime corresponds in our fits to a linear evo- 
lution of the mode-coupling parameters with temperature. Approaching the 
MCT critical transition, the effective vertices depart from the linear behavior 
and exhibit an unexpected behavior : they follow the transition line with- 
out crossing it. The avoidance of the transition has been known for a long 
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time and is generally associated with the occurrence of so-called activated 
processes that are not included in the theory. However, in our picture, these 
competitive processes seem to be efficient well above the critical tempera- 
ture predicted within the MCT formalism, in a domain where activated or 
hopping processes are not supposed to be important. An immediate conse- 
quence of this is that the universal, critical predictions of the theory should 
be unobservable experimentally. 

That the MCT of the glass transition needs be completed is not a new 
idea, but it is a priori a difficult task. If the picture given here is valid, one 
would expect from a successful extension of MCT that it restores the quasi- 
linear variation of the effective vertices by introducing new ones that become 
important only at low temperature. In this respect, note that our results for 
the critical parameters extrapolated from high temperatures (where the nec- 
essary corrections to ideal MCT are thought to be of minor importance) do 
not allow to recover those from the so-called "extended mode-coupling the- 
ory" that has been proposed as an improvement of the ideal mode-coupling 
theory and was tested in its scaling form on salol and CKN.i 

At last, it is necessary to point some remaining difficulties, that are of 
great interest but are out of the scope of this study. Of critical importance 
is the nature of the relaxation probed by the depolarized light scattering, 
density fluctuations or molecular reorientations.il Fortunately, we find that 
the choice for the expression of the scattered intensity as a function of the 
schematic correlators has no influence on the results of fit. But it asks the 
question of the dynamics of molecular orientational processes in the super- 
cooled liquid phase and its role in the structural relaxation of molecular 
systems. Some evidence exists in simulations that this dynamics is the same 
as for density fluctuations and that the two processes are coupled.il The 
very good results of the schematic fits are also in favor of such a conclusion. 
Nevertheless, this remains to be demonstrated theoretically. Thus, in order 
to test the relevance of this approach from the experimental point of view, 
analysis of data from different spectroscopic techniques are necessary. With 
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this goal, coherent neutron scattering spectra (that probe directly the dy- 
namic structure factor) and far infrared absorption data (that probe dipolar 
reorientations, thus molecular rotations) are now under study 
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Note 

After submission of this work, a paper was sent to us by W. Gotze, which 
reports schematic calculations on DLS spectra of glycerol.il The model used 
in this paper and the results obtained will be discussed and compared to ours 
in further publication.il 



23 



References 

[1] C.A. Angell, p. 3 in Relaxation in Complex Systems, edited by K.L. Ngai 
and G.B. Wright (Office of Naval Research, Washington, 1984); C.A. 
Angell, J. Phys. Chem. Solids 49, 863 (1988); C.A. Angell, in Proceedings 
of the first Tohwa University International Symposium, Fukuoka, Japan, 
1991, AIP Conf. Proc. 256, edited by K. Kawasaki, T. Kawakatsu and 
M. Tokuyama (AIP, New York, 1992). 

[2] W. Gotze, Philos. Mag. 43, 219 (1981) ; E. Leutheusser, Phys. Rev. 
A 29, 2765 (1984); W. Gotze, p.287 in Liquids, Freezing and Glass 
Transition (Les Houches, 1989) , edited by J. P. Hansen, D. Levesque 
and J. Zinn- Justin (North-Holland, Amsterdam, 1991) and references 
therein; W. Gotze and L. Sjogren, Rep. Prog. Phys. 55, 241 (1992). 

[3] T. Franosch, M. Fuchs, W. Gotze, M.R. Mayr and A.P. Singh, Phys. 
Rev. E, in print. 

[4] C. Alba-Simionesco and M. Krauzman, J. Chem. Phys. 102, 15 (1995). 

[5] C. Alba-Simionesco, V. Krakoviack, M. Krauzman, P. Migliardo and F. 
Romain, J. Raman Spectroscopy 27, 715 (1996). 

[6] V. Krakoviack, C. Alba-Simionesco and M. Krauzman, p. 373 in Non- 
equilibrium phenomena in supercooled fluids, glasses and amorphous ma- 
terials, Pisa, Italy, 1995, edited by M. Giordano, D. Leporini and M.P. 
Tosi (World Scientific, Singapore, 1996). 

[7] G. Li, W.M. Du, A. Sakai and H.Z. Cummins, Phys. Rev. A 46, 3343 
(1992). 

[8] H.Z. Cummins, W.M. Du, M. Fuchs, W. Gotze, S. Hildebrand, A. Latz, 
G. Li and N.J. Tao, Phys. Rev. E 47, 4223 (1993). 

[9] J. Toulouse, G. Coddens and R. Pattnaik, Physica A 201, 305 (1993). 



24 



[10] C. Dreyfus, M.J. Lebon, H.Z. Cummins, J. Toulouse, B. Bonello and 
R.M. Pick, Phys. Rev. Lett. 69, 3666 (1992). 

[11] Y.W. Yang and K.A. Nelson, Phys. Rev. Lett. 74, 4883 (1995); Y.W. 
Yang and K.A. Nelson, J. Chem. Phys. 103, 7732 (1995). 

[12] G. Li, W.M. Du, X.K. Chen, H.Z. Cummins and N.J. Tao, Phys. Rev. 
A 45, 3867 (1992). 

[13] F. Mezei, W. Knaak and B. Farago, Phys. Rev. Lett. 58, 571 (1987); 
W. Knaak, F. Mezei and B. Farago, Europhys. Lett. 7, 529 (1988); F. 
Mezei, J. Non-Cryst. Solids 131-133, 317 (1991). 

[14] G. Li, W.M. Du, J. Hernandez and H.Z. Cummins, Phys. Rev. E 48, 
1192 (1993). 

[15] Y.W. Yang and K.A. Nelson, J. Chem. Phys. 104, 5429 (1996). 

[16] P. Lunkenheimer, A. Pimenov and A. Loidl, submitted to Phys. Rev. 
Lett., 10/96. 

[17] W. Steffen, A. Patkowski, G. Meier and EW. Fisher, J. Chem. Phys. 
96, 4171 (1992); W. Steffen, G. Meier, A. Patkowski and EW. Fisher, 
Physica A 201, 300 (1993); W. Steffen, A. Patkowski, H. Glaser, G. 
Meier and EW. Fisher, Phys. Rev. E 49, 2992 (1994). 

[18] A. Patkowski, W. Steffen, G. Meier and EW. Fischer, J. of Non-cryst. 
Solids 172-174, 52 (1994); W. Steffen, B. Zimmer, A. Patkowski, G. 
Meier and EW. Fischer, J. of Non-cryst. Solids 172-174, 37 (1994). 

[19] J. Jackie, p. 135 in Amorphous solids : low-temperature properties, edited 
by W.A. Phillips (Springer, Berlin, 1981). 

[20] V. Malinovski and A.P Sokolov, Solid State Commun. 57, 757 (1986); 
A. P. Sokolov, in Dynamics of disordered materials II, edited by A.J. 
Dianoux, W. Petry and D. Richter (North- Holland, 1993); U. Buchenau, 



25 



Europhys. News 24, 77 (1993); V.N. Novikov, E. Duval, A. Kisliuk and 
A.P. Sokolov, J. Chem. Phys. 102, 4691 (1995). 

[21] E. Duval, N. Garcia, A. Boukenter and J. Serughetti, J. Chem. Phys. 99, 
2040 (1993); T. Achibat, A. Boukenter and E. Duval, J. Chem. Phys. 
99, 2046 (1993). 

[22] H.Z. Cummins, G. Li, W.M. Du, R.M. Pick and C. Dreyfus, Phys. Rev. 
E 53, 896 (1996); erratum in Phys. Rev. E 55, 1232 (1996). 

[23] S.P Das and G.F. Mazenko, Phys. Rev. A 34, 2265 (1986); B. Kim, 
Phys. Rev. A 46, 1992 (1992); J. Yeo and G.F. Mazenko, Phys. Rev. E 
51, 5752 (1995). 

[24] T. Geszti, J. Phys. C 16, 5805 (1983) 

[25] L. Sjogren, Phys. Rev. A 33, 1254 (1986). 

[26] J. Bosse and U. Krieger, J. Phys. C : Solid State Phys. 19, L609 (1986); 
U. Krieger and J. Bosse, Phys. Rev. Lett. 59, 1601 (1987); U. Krieger 
and J. Bosse, Z. Phys. Chem. 156, 97 (1988). 

[27] N.J. Tao, G. Li, X. Chen, W.M. Du and H.Z. Cummins, Phys. Rev. A 
44, 6665 (1991). 

[28] M.J. Stephen, Phys. Rev. 187, 279 (1969). 

[29] V. Krakoviack, C. Alba-Simionesco, in preparation. 

[30] M. Fuchs, W. Gotze, I. Hofacker and A. Latz, J. Phys. : Condens. 
Matter 3, 5047 (1991). 

[31] P.J. Davis and P. Rabinowitz, in Methods of Numerical Integration (Aca- 
demic press, New York, 1975). 

[32] W.H. Press, S.A. Teukolsky, W.T. Vettering and B.P Flannery, in Nu- 
merical Recipes (Cambridge university press, New York, 1985). 



26 



[33] X.C. Zeng, D. Kivelson and G. Tarjus, Phys. Rev. E 50, 1711 (1994); 
P.K. Dixon, N. Menon and S.R. Nagel, Phys. rev. E 50, 1717 (1994); 
H.Z. Cummins and G. Li, Phys. Rev. E 50, 1720 (1994). 

[34] W. Kob and H.C. Andersen, Phys. Rev. Lett. 73, 1376 (1994); W. Kob 
and H.C. Andersen, Phys. Rev. E 51, 4626 (1995); W. Kob and H.C. 
Andersen, Phys. Rev. E 52, 4134 (1995); M. Nauroth and W. Kob, 
Phys. Rev. E 55, 657 (1997). 

[35] B. Bernu, J. P. Hansen, Y. Hiwatari and G. Pastore, Phys. Rev. A 36, 
4891 (1987); J.N. Roux, J.L. Barrat and J. P. Hansen, J. Phys. : Con- 
dens. Matter 1, 7171 (1989); J.L. Barrat and A. Latz, J. Phys. : Con- 
dens. Matter 2, 4289 (1990). 

[36] W. G6tze and L. Sjogren, J. Phys. C : Solid State Phys. 21, 3407 (1988). 

[37] W. Gotze and L. Sjogren, J. Phys. : Condens. Matter 1, 4183 (1989). 

[38] G.F. Signorini, J.L. Barrat and M.L. Klein, J. Chem. Phys. 92, 1294 
(1990). 

[39] T. Franosch, W. Gotze, M.R. Mayr and A.P Singh, Phys. Rev. E 55, 
3183 (1997). 



27 



Table 1: Parameters for MCT fits to salol data (T m = 315 K, T g = 218 K). 



technique 


T C (K) 


A 


references 


DLS|, ideal MCT 


256 ± 5 


0.70 





DLS a , extended MCT 


250 


0.73 


1 


neutron scattering 


263 ± 7 




1 


Brillouin scattering 


275 ± 10 







viscosity 


265 ± 5 






ISLS0 


266 ± 1 








"depolarized light scattering 
impulsive stimulated light scattering 



Table 2: Parameters for MCT fits to CKN data (T m = 433 K, T g = 333 K). 



technique 


T C (K) 


A 


references 


DLS0, ideal MCT 


378 ± 5 


0.81 




12 




DLS a , extended MCT 


378 


0.85 




1 


neutron scattering 


368 ± 5 






13 


Brillouin scattering 


375 ± 5 






14 




ISLSQ 


378 ± 2 






15 




Dielectric spectroscopy 


375 


0.76 




16 



"depolarized light scattering 
impulsive stimulated light scattering 



28 



Table 3: Values of the temperature-independent parameters of the schematic 
model 





salol 


CKN 


fio/2vr (GHz) 


1230 


2000 


r 


0.11 


0.3 


^i/2tt (GHz) 


160 


900 


Ti 


5.2 


2.6 


r 


37 


10 


7 


18.6 


2.4 



Table 4: Parameters from the schematic-model fits to salol and CKN data 





salol 


CKN 


T e (K) 


257 ± 5 


388 ± 5 


Xht 


0.69 ± 0.02 


0.73 ± 0.02 


Xvert 


0.61 ± 0.02 


0.69 ± 0.02 
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FIGURE CAPTIONS 



Figure 1: Experimental depolarized light scattering spectra (for readabil- 
ity, some spectra are omitted and the data are reported every 2 points) 
and schematic fits, a) salol (the last two curves at 233 K and 253 K 
are obtained with a 90 GHz low frequency cut-off; the spectra at 243 
K, 273 K, 293 K, 313 K and 333 K are omitted), b) CKN (the spectra 
at 383 K, 403 K, 423 K, 443 K and 468 K are omitted). 

Figure 2: Individual contributions to the fitted DLS spectra of each of 
the two correlators of the schematic model in a semilogarithmic plot, 
a) salol at 333 K. b) CKN at 433 K. Note that the Xo" contribution is 
negligible below 50 GHz for salol, but not for CKN. This explains the 
equivalence between a scaling-equation fit and our schematic calcula- 
tions in the case of salol. 

Figure 3: Vertices V\ and i>2 of the schematic model as functions of tem- 
perature for salol and CKN. The straight lines correspond to the linear 
extrapolation described in the text. 

Figure 4: Two-dimensional (t>i,i>2) parameter space of the schematic mo- 
del and trajectories followed by salol and CKN with temperature. 

Figure 5: Effective and local exponent parameters (resp. X e ff and A^ oc ) as 
functions of e (see text for definition) for the one-correlator F 12 -model. 

Figure 6: Effective temperature-dependent, effective temperature-indep- 
endent, and local exponent parameters (resp. A e //, Xgiobai, and A; oc ) as 
functions of temperature for supercooled salol above T c . 

Figure 7: Individual contributions to the fitted DLS spectra of each of 
the two correlators of the schematic model at the closest temperatures 
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higher than T c (the scaling factor is such that the minima of suscep- 
tibilities are equal), a) salol at 263 K. b) CKN at 393 K. Note that 
the minima of both contributions (when they exist) occur at differ- 
ent frequencies, in disagreement with the asymptotic predictions of the 
mode-coupling theory. 
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